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Abstract 

Effective potentials of the relativistic ma 6 order correction for the 
ground state of the Coulomb two-center problem are calculated. They 
can be used to evaluate the relativistic contribution of that order to the 
energies of hydrogen molecular ions or metastable states of the antipro- 
tonic helium atom, where precision spectroscopic data are available. In 
our studies we use the variational expansion based on randomly chosen 
exponents that permits to achieve high numerical accuracy. 

1 Introduction 

In recent years several experiments on precision spectroscopy of three-body light 
atomic and molecular systems become available [TJ [2j E] . That is a great chal- 
lenge to theorists, since the quantum electrodynamics for the few-body bound 
state problem is not that well elaborated as for atoms comprised of two particles. 

On the other hand, numerical variational solution for the nonrelativistic 
Coulomb three-body problem have reached very high precision. The ground 
state of the molecular ion has been calculated with as much as 30 digits [4 
(for other examples of nonrelativistic calculations see references therein). 

Accurate knowledge of both theoretical and experimental transition frequen- 
cies has the potential of improving m e /m p , m e /mp, etc mass ratios [5J HI [7J- 
In order to achieve this goal the important step is evaluating the relativistic 
and radiative corrections for the binding energies. This can be systematically 
performed using series expansion of the binding energy in terms of the coupling 
constant, a. The leading order relativistic and radiative corrections (RooCt 2 , 
Rood 3 ) have been obtained in [SJ [5] . The main aim of this work is to calculate 
a contribution of order ma e (or i? oc ,a ) due to relativistic corrections for the 
bound electron in the nonrecoil limit. The radiative corrections to this order are 
known analytically jlOj . A sum of these two contributions constitutes a com- 
plete ma 6 order correction for a one electron three-body system in a nonrecoil 
limit. 



1 



2 Nonrelativistic Hamiltonian and Variational 
Expansion 

In what follows we consider the Coulomb two-center problem with the nonrel- 
ativistic Hamiltonian 



2m 



r\ r 2 



(1) 



where r\ and r 2 are the distances from an electron to nuclei 1 and 2, respectively. 
Atomic units, m e = e = h=l, are used throughout. 

In order to get a precise solution for the Schrodinger equation 



2m 



V 



*o(r) = ^o*o(r), 



(2) 



we use the variational approach. A variational expansion for the wave function 
of the ground state of electron is taken in the form (Z\ ^ Z 2 ): 



*(r) 



e imip r \m\ 



y * (J. e -a i r 1 -/3ir 2 
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where r is a distance from an electron to the z-axis and 
1 



2R 



2r\r\ + 2r\R? + 2r\R? -r\-r\- R 4 . 



For Z\ = Z 2 the variational wave function should be symmetrized 

00 

'y ' (J,^ e —ctiri—0tr2 _|_ g— ftri— a*r2 ^ 



(4) 



i=i 



where (+) is used to get a gerade electronic state and (— ) is for an ungerage state, 
respectively. Parameters an and pi are generated in a quasi-random manner 



(A 2 -A 1 ) + A 1 



(5) 



Here \_x\ designates the fractional part of x, p a is a prime number, an interval 
L4i, A2] is a real variational interval, which has to be optimized. Parameters (3i 
are obtained in a similar way. 

The Schrodinger equation expressed in coordinates, r%, r 2l <p, has a form: 



T\T 2 
Z X Z 2 



d d - — 



n 



(6) 
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In order to get accurate results we use several sets of basic functions of Eq. 
(or ((4])) (in a spirit of [11]) with different values of variational parameters: three 
sets for small values of internuclear distance R, and two sets for intermediate 
and large values of R, respectively. Total number of basis functions varies 
from N — 110 to N — 250. Using the variational method described above 
the nonrelativistic energy has been obtained with accuracy of more than 20 
significant digits. In these calculations arithmetics of sextuple precision (about 
48 decimal digits) has been used. Earlier, in a previous paper [T3] using the same 
variational approach the effective potentials of the Breit-Pauli Hamiltonian have 
been calculated with up to 10 significant digits. 

The matrix elements, which appear in this and following sections, have been 
evaluated analytically, for further details see Appendix A. 



3 Breit— Pauli Hamiltonian 

The leading order relativistic correction of order ma can be expressed by the 
Breit-Pauli Hamiltonian: 



P 4 1 

H B = + Zi47rf(ri) + Z 2 ^5{r 2 )Y 

8m J 8m z 



[n x p] [r 2 x p] ' 

Z l~^ 2 3~~ + Z2 ~^ TT~ S ■ 



(7) 



The spin-dependent term (the last term in Eq. of the Breit-Pauli Hamil- 
tonian H S B ° does not contribute to the energy at this order, but should be taken 



4 Relativistic corrections of ma order 

In this section we will assume for simplicity of notation that Hb = H^, the 
scalar term of Hb- The spin-orbit term will be treated separately. 

Our method is based on the effective Hamiltonian approach, which is closely 
related to ideas of [14l[T5], where a case of atomic hydrogen has been considered. 



4.1 Formal expressions 

The energy displacement due to the ma 6 order relativistic corrections can be 
expressed as: 

A£< 6 ) = (H B Q(E - H^QHb) + (# (6) ) . (8) 
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Here Q = I— \ijjo)(ipo\ is a projection operator. The effective Hamiltonian 
corresponding to this contribution has a form [14j 



16m 5 8m 3 16m 4 

5 



(9) 



128m 4 v ; 64m 4 

where Si = ~Z i r l /rf and pi = ZiSfa) (AV = 4np). 
Separately, both contributions of (|8|) are divergent. 

4.2 Removing divergences from the second order contri- 
bution 

The second order perturbation term 

AE B 2} = (H B Q(E - H )- X QH B ) (10) 
can be evaluated by obtaining the wave function $b as a solution of equation 

(Eq - H )^b = {H B - {H B ))* , (11) 

_____ r n| 

then the correction to the energy from Eq. (jTUJ) can be expressed as AE B = 
(9 \(H B -(H B ))\9 B ). 

It is known that a formal second order expression of the type 

(^o\S 3 (ri)Q(E - flo) _1 Q<5 3 (ri)|*o) 
is divergent. Similarly p 4 |^o) behaves as 

-2mp 2 V\* ) ~ 8TTm[Z 1 5(r 1 ) + Z 2 5{r 2 )}\^ Q ) 

at small values of ri , r 2 ■ 

In order to understand how to overcome this problem, let us consider first 
the atomic hydrogen S-state case. Let ^ B be a solution of equation 



En- 



f__Z_ 

2m r 



*<f> = Q 



p Zir 



8m 3 2m 2 



8(t) 



4 H \ (12) 



separating the dominant behaviour at small r, the solution of above equation 
can be presented as 



Amr 



where \P B is a less singular function, "J^ ; ~ lnr at r — > 0. 

Coming back to the two center problem, let us try to separate the singular 
part of the Breit-Pauli wave function solution $b as follows: 

*b = C/*o + *b. (13) 
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where U = —^V. Substituting of Eq. (fT5|) into Eq. (JTUJ) modifies our equations, 
AE { B 2} = (tfoK-ff* - (H B ))U\*o) + (9 \(H B - (H B ))\* B ) , (14a) 

and $b is a solution of 

(E - H )* B = -{Eo - H )Uy + (H B - (H B ))9 . (14b) 

However, the last term in Eq. (|14ap is still divergent, the singularity from the 
left-hand side of this term should be eliminated as well: 



*o\(H B -(H B ))\V B/ 

= (VoKHB-iHB^iEo-Ho)- 1 

x[-(E -H )U + (H B -(H B ))}\^ ) (15) 
= (* | -U(E Q -H )U + U(H B -(H B ))\y ) 

* fl | - (E -H )U+ {H B -(H B ))\9 ) 



These transformations as can be seen are equivalent to the one used by Pachucki 
in Pi]: 

'H' B =H B ~ (E Q -H )U~U(E Q -H Q ) 

(HbQ(E -H )- x QHb) = {H' B Q{E -H )- l QH' B ) (16) 

+ (UH B + H B U) - 2 (U) (H B ) - (U(Eo-H )U) 

The last three terms can be recast in a form of a new interaction modifying the 
Hamiltonian H^: 

H'^ = (U H B + H B U) - 2U(H B ) - U(E a - H )U, 

(H' i6) ) = (P 4 V+V P *) - ( V[ P1+P2 ] > (17) 

A respective transformation of the Breit-Pauli operator is 

H' B = -^3 + ^[ZiSfr) + Z 2 6(r 2 )] 

1 ;(£ 1 +£ 2 )V + 2U(H -E a ) 



(18) 



4m 2 

Finally, the energy shift can be written 

AE^ = (H' B Q(E - H )- l QH' B ) + (h 1 ^) + (h^) (19) 



And now divergent terms are gathered together into the modified effective 
Hamiltonian H ^ + H'^ . 
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4.3 Removing divergences from the modified effective Ha- 
miltonian 




G 



5 Results and Conclusion 



Equations (fl9|) . (|2Tj) . and (|23| have been evaluated numerically using the vari- 
ational expansion ©-JI]). 

In Tables 1 and 2 the relativistic corrections of order ma 6 for the ground state 
of electron for symmetric Z\ = Z% = 1 (H^ molecular ion) and asymmetric (the 
antiprotonic helium atom) Z\=2 and Z^ = — 1 cases are presented as functions 
of a bond length. The accuracy of obtained results for AE 1 * 6 * 1 is estimated as all 
digits indicated in the Tables. Figures 1 and 2 shows the "effective" potentials 
for these two cases. Dashed lines are the radial wave function for the ground 
and first vibrational S-states of and the (36, 34) state of the 4 He + p atom. 

The last table gives a comparison of our calculations with the earlier ones, 
which demonstrate superiority of the newly obtained results. However, we want 
to draw attention to the approach of [18] , which, to our opinion, is very promis- 
ing and allows to evaluate the ma 8 order corrections as well using the same first 
order perturbation wave function with a rather high accuracy. Presumably, the 
not very high precision of this particular calculation is connected with the Gaus- 
sian basis set, which is not very suitable for description of the zero-order and 
first-order solutions. 
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A Analytical evaluation of the matrix elements 

The calculation of the matrix elements is reduced to evaluation of integrals of 
the type 

Tim (a, /3)= [ r[- l r™- l e-^-^d?T. (24) 



Integers {l,m) are, in general, non- negative, but in case of singular matrix ele- 
ments one of the indices can be negative. 
The function Tqq can be easily obtained 

Att e~P R — e~ aR 

Too(a,0,R) = - a2 _ p2 , (25) 

where R is the distance between nuclei, then Ti m (a, f3; R) for non-negative (I, m) 
may be generated from (f25"|) by means of relation 

d V / d 



r lm (a,0;R)= (-—) I-—) r 00 (a,{3,R). (26) 
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Integral r_i i o(o;, /?; R) is expressed by 

T_ h0 (a,p;R) = — le^ExiRia + /3)) + e""* In R(a + p) 

R/3 I (27) 

-e^[Ei(i?(a - /?)) + In iJ(a -/?)]}• 

Worthy to note that a function in square brackets is analytic when argument is 
zero. Integrals r_i j7n are generated from r_i.o similar to 



~dBJ T -w( a >P> R )- ( 28 ) 
Function Ei(z) encountered in (f27|) is the exponential integral function [20] : 



/•OO 

Ei(z) = r(0,z) = / t^e^dt. 

J z 



B Relations between divergent matrix elements 

In this section we will assume that V is regularized in some or other way and 
p = pi + p 2 (AV — A-np) is a smooth function of space variables, a distribution 
of charge in space. The left- and right-hand side functions inside the brackets 
represent the same wave function \l/o, the solution of the Schrodinger equation 
((2). Then by using commutation relations and integration by parts one gets 

(Vp 2 V) = (V 2 p 2 ) - Air (V( Pl +p 2 )) + 2(V(£ 1 +£ 2 )V) , (29a) 
(Vp 2 V) = ((£ 1 +£ 2 ) 2 ) - 2 {V(£ 1 +£ 2 )V) + ( P V 2 p) , (29b) 

4ir(V( Pl +p 2 )) = - ((£!+£ 2 ) 2 ) + 2 (V(£ 1 +£ 2 )V) . (29c) 
Subtracting Eq. (|29c[) from Eq. (|29a|) . we have 

(Vp 2 V)=({£ l+ £ 2 ) 2 ) + (V 2 P 2 ), (30a) 

then summing up (|29b[) and (|29c|) and taking into account (|30a|) : 

4ir(V( Pl +p 2 )) = - ((£ 1 +£ 2 ) 2 ) - (V 2 p 2 ) + ( P V 2 p) . (30b) 

From the last two formulas with a systematic use of 

p 2 ^o = 2m(E -V)^o, 

the final expressions for the divergent terms in Eq. (|2"D|) may be obtained: 

f (p 6 ) _{(£ 1 + £ 2 ) 2 ) (V 3 ) t 3E (V 2 ) 3E 2 (V) , E% 



16m 5 4m 3 2m 2 2m 2 2m 2 2m 2 ' 

(p^V+Vp 4 ) = ((£i+£ 2 ) 2 ) + CO _ E o (V 2 ) + Eq (V) 
128m 4 32m 3 16m 2 8m 2 16m 2 ' 

*(v(pi+P2)) = ((£i+£ 2 ) 2 ) (v 3 ) e (v 2 ) + ( P v 2 p) 

2m 3 8m 3 4m 2 4m 2 8m 3 



(31) 
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After summing up all the terms in (|20p the regularization can be removed. 
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Figure 1: Adiabatic "effective potential" for the ma 6 order relativistic correc- 
tion, AE^ (Zi =Z 2 = 1) and radial wave functions of the first two vibrational 
states (L = 0) of the molecular ion. 




Figure 2: Adiabatic "effective potential" for the ma 6 order relativistic correc- 
tion, AE^ (Zi = 2, Z 2 = -1) and a radial wave function of the (L = 34, v = 2)) 
state of the 4 He + p atom. 
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Table 1: Adiabatic potentials of the nonrelativistic energy and the ma 6 order 
relativistic correction, AE^\ for Z\ = Z2 = 1. 
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Table 2: Adiabatic potentials of the nonrelativistic energy and the ma 6 order 
relativistic correction, AE^ 6 \ for Z\ — 2, Z2 = — 1. 
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Table 3: Comparison with earlier calculations at a bond length R = 2.0 [Bohr]. 
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